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ABSTRACT 

Taking into account the mixing effects between left- and right-handed top-squarks, 
we calculate the genuine supersymmetric eletroweak correction to top quark pro- 
duction at the Tevatron in the minimal supersymmetric model. The analytic ex- 
pressions of the corrections to both the parton level cross section and the total 
hadronic cross section are presented. Some numerical examples are also given to 
show the size of the corrections. 
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1. Introduction 



The top quark has been discovered by CDF and DO collaborations at the Teva- 
tron [1]. The mass and production cross section are found to be 176 ± S(stat) ± 
10(syst) (l99±H(stat) ± 22(syst)) GeV and 6.8±H (6.4 ± 2.2) Pb by CDF ( DO ) 
collaboration. The comparison of the theoretical calculation of the top quark produc- 
tion cross section with experimental results is necessary in the test of the mechanism 
by which top quarks are produced. Within the framework of the Standard Model 
(SM) the next-to-leading-order calculation for the QCD processes was completed sev- 
eral years ago[2]. Recent works [3] have extended those results with the inclusion of 
the exact order corrected cross section and the resummation of the leading soft 
gluon corrections in all orders of perturbation theory. The cross section was predicted 
to be a t t(m t = UQGeV) = 4.79t[J;gP6[3]. The latest results given by Berger [4] was 
Vtt{™>t = l75GeV) = 5.52t%%Pb. The one-loop electroweak corrections to the cross 
section were found to be only a few percent [5]. Therefore, the results of theoretical 
prediction in the SM are almost consistent with the experimental results within the 
region of the errors. 

Since the corrections to top quark production cross section above 20% are poten- 
tially observable at the Tevatron, it is tempting to calculate the radiative corrections 
arising from the new physics beyond the SM. In the minimal supersymmetric model 
(MSSM), the Yukawa correction and supersymmetric QCD correction were calculated 
in Refs.[6,7]. These corrections cannot reach the observable level for experimentaly 
allowed parameter values. The genuine supersymmetric electroweak corrections of 
order amj/m^, which arise from loops of chargino, neutralino and squark, have also 
been calculated by us in Ref. [8] and its erratum [9], where we neglected the mix- 
ings between left- and right-handed squarks and assumed the mass degeneracy for all 
squarks. In such a simple case, the analytic results were quite simple and the numeri- 
cal size of the corrections could not reach the observable level for squark mass heavier 
than 100 GeV. However, due to the possible significant mixing effects between left- 
and right-handed top-squarks, which is suggested by low-energy supergravity models 
but is completely general [10], the mass splitting between the two mass eigenstates 
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of top squark may be quite large. The supersymmetric electroweak corrections may 
be sensitive to top-squark mixing effects. 

In this paper, taking into account the mixing effects between left- and right- 
handed top-squarks, we present the genuine supersymmetric eletroweak correction to 
top quark production at the Tevatron in the minimal supersymmetric model. In Sec. 2, 
we briefly overview top-squark mixing. In Sec. 3, we give the analytic expressions of the 
corrections to both the parton level cross section and the total hadronic cross section. 
In Sec. 4 we present some numerical examples to show the size of the corrections. 

2. Top-squark mixing 

The mass matrix of top-squarks takes the form [10] 

~ Lm = (tl ( m t M L R ) 0* ) 

ml = Mf T +m 2 t + cos(2/3)(^ - - sin 2 9 W )M 2 Z 



m 2 



= M 2 +m 2 + - cos(2/3) sin 2 6 W M 2 Z 



tR t R 1 * 3 

M LR = ficot (3 + A t M (I) 

where M~ , M? are the soft SUSY-breaking mass terms for left- and right-handed 
top-squarks, \i is the coefficient of the H\H 2 mixing term in the superpotential, A t M 
is the coefficient of the dimension-three trilinear soft SUSY-breaking term t L t R H 2 , 
and tan/5 = v 2 /vi is the ratio of the vacuum expectation values of the two Higgs 
doublets. 

The mass eigenstates of top-squark are obtained by 

( h \ _ ( t L \ _ ( cose sine \ ( t L \ 

\t 2 )~ \t R )-{ -sine cose ){t R ) [Z) 

and the masses of t± t2 are given by 

R Ml RT=( < ° ) (3) 



*2 



The expressions of 6 and are given by 



tan 20 = (4) 



m L = 2 



m~ — ra- 
th tR 

1 



m h + m in =F V ( m l - m V 2 + 4M LR m t (5) 



For sbottoms, since we neglect the mixing between left- and right-handed sbot- 
toms, we have 

m L = m t L ,n =m 2 b + Ml R ±co S (2P)(Tl R -Q bS m 2 9 w )M 2 z (6) 

where T^ R — — |, and Qb = — |. M~ hh are the soft SUSY- breaking mass terms for 
left- and right-handed sbottoms. 

3. Analytical expression of the correction 

At the Tevatron, the top quark is dominantly produced via quark-antiquark anni- 
hilation [11]. The genuine supersymmetric electroweak correction of order am^/rayy 
to the amplitude is contained in the correction to the vertex of top-quark color cur- 
rent. The relavent Feynman diagrams are shown in Fig.l in Ref. [7]. The Feynman 
rules can be found in Ref. [12]. In our calculation, we use dimensional regulariza- 
tion to regulate all the ultraviolet divergences in the virtual loop corrections and we 
adopt the on- mass-shell renormalization scheme [13]. The renormalized amplitude for 
qq — > ti can be written as 

M ren = M + 5M, (7) 

where M is the amplitude at tree-level and SM is the correction to the amplitude, 
which are given by 

M = v{p 2 ){-ig s T A Y)u{ Pl ) u(p 3 )(-ig s T A Y)v(p4), (8) 

s 

5M = v(p 2 )(-ig s T A Y)u( Pl ) u(p 3 )5A»v(p 4 ). (9) 

s 

Here, pi,p 2 denote the momenta of the incoming partons, and P3,p^ are used for 
outgoing t quark and its antiparticle. s is center-of-mass energy of parton level 
process. 5A^ stand for the genuine supersymmetric electroweak correction to the 
vertex of top-quark color current, which are given by 

+k»F 3 + k^ l5 F 4 + ik v a^F b + tk u a^ l5 F 6 ] , (10) 
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where = |[7 M , 7"] and the form factor Fj are obtained by 



F = F c + F n 



(11) 



where Ff and F" arise from chargino and neutralino digrams, respectively. Ff are 
given as 



Ft = 



E V ^ V f2 c 24 + m 2 t {c n + c 21 ) + £-B x + m 2 B[)(m t , M,-, m~ b ) 
3=1,2 



n 



n 



c 24 + -Biimt.M^mi) 



E ^ 

J=l,2 

-m* E ^2^/2 (C21 - 2c 23 ) , 



i=i,2 



o m « E ^2^/2 (C21 + 4c 22 - 4c 23 ) , 



.7=1,2 



E ^2^/2 (en + C21) , 



j=l,2 



-m t E ^2^/2 (en + c 2 i - 2c i2 - 2c 23 ) , 



j=l,2 



(12) 

(13) 

(14) 

(15) 

(16) 

(17) 
(18) 



where the functions Cy(— p 3 ,p 3 + P4, M,-, m^, m^) and I?i are the 3-point and 2-point 
Feynman integrals [14]. The chargino masses Mj and matrix elements Vy depend on 
parameters M, //,tan/3, whose expressions can be found in Ref. [12]. _Bq ;1 are defined 
by 

dB 0>l (p,m 1 ,m 2 ) 



dp 2 



Fp are obtained by 



F? = F- 1 + F- 2 + F t s (fori = 1,2) 
if\» = F^+Ft (fori = 3, 4, 5, 6) 



(19) 

(20) 
(21) 



F{ and F| are given by 



j'=i 



F i S = E o^'4iV; 4 5i(m () M 0j , m fl ) + Si(mt, M 0j , m f J 



+m 2 N j4 N* 4 [B[(m t , M 0j ,m h ) + ^(m t , M 0j , m t ~ 2 ) 
+m t M 0j N j4 N j4 sin(2#) [^(m t , M 0j , m t ~ 2 ) - B' (m t , M 0j , m- tl )] } (22) 
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F 2 = E 0^4^*4 cos(2^) [Si(mt, M 0j , m h ) - B x {m t , M 0j , m h ) 



(23) 



j'=i 

F/ 1 are given by 



F* 1 = 



C 2 4 + ^t(Cll + C 2 lj 



if 



H 1 



Fl 1 



sm(20)N j4 N j4 m t M Oj (co + c u )} (24) 

(25) 
(26) 



E {^4^4 
J'=l 
4 

E N j4 N* 4 c 24 cos(2#) 
j'=i 

4 ri 1 

2 -m t N j4 N* 4 (c 21 - 2c 23 ) + - sin(2^)A^-4^4M 0j (2c 12 - c n ) 
1 4 

- cos(2fl)m t ]T A^. 4 Ag 4 ( C21 + 4c 22 - 4c 23 ) (27) 

4 r 1 1 1 

J2 -^rn t N ]A N* A ( Cll + c 21 ) + - sin^^A^c^Co + c n ) (28) 
j=i L ^ I 

-- cos(2#)m t £ N j4 N* 4 (c u - 2ci 2 + c 21 - 2c 23 ) (29) 



where P3,P3 + P4, M j, m^, m^) are the 3-point Feynman integrals [14]. The 
neutralino masses Mqj and matrix elements Nij are obtained by diagolis the matrix 
Y [12]. Giving the values for the parameters M, M', /x, tan/3, the matrix iV and 
AT D can be obtained numerically. Here, the parameters M, M' are the masses of 
gauginos corresponding to SU(2) and U(l), respectively. With the grand unification 
assumption, i.e. SU(2) x U(l) is embedded in a grand unified theory, we have the 



relation M' = \^M. F- 2 are given by 



T?t 2 _ pt 



sin(20)-+- sin( 2 0),cos( 2 0)^- cos(20),m£ 



*2 



(30) 



The renormalized cross-section for parton level process qq — > tt are given by 

a(s) = ct° + Act, 



(31) 



with 



27l2 



A(s + 2m t 2 ), 



Act 



l 0i 



g 2 m 2 t 



9s 3 ^ 327r 2 m^ sin 2 (3 



-F ± s(s + 2m 2 ) + 2F 5 m t s 2 



(32) 
(33) 



where /5 t — y 1 — Amf/s. 

The hadronic cross section is obtained by convoluting the subprocess cross section 
&ij of partons i,j with parton distribution functions ff(xi, Q), ff(%2, Q), which is 
given by 



with 



dLij r 1 dxi 



f -^\jf(x 1 ,Q)ff(r/x 1 ,Q) + (A^B)} (35) 

Jr Xi l j j 



rfr Jt x 1 

In the above the sum runs over all incoming partons carrying a fraction of the proton 
and antiproton momenta (p lj2 = ^1,2-^1,2)5 v^S* = 1.8 TeV is the center-of-mass energy 
of Tevatron, r = x\Xi and r = Amf/S. As in Ref.[3], we do not distinguish the 
factorisation scale Q and the renormalisation scale ji and take both as the top quark 
mass. In order to compare our results with the Yukawa corrections in Ref. [6], we use 
the same parton distribution function as in ref. [6], i.e. the Morfin-Tung leading-order 
parton distribution function [15]. 

4. Numerical examples and discussion 

In the numerical examples presented in Figs. 1-3, we fixed M = 200GeV, fi = 
-lOOGeV and used the relation M' = f^M to fix M' . Also we assumed M- tR = 
Mf L = M~ bL which depend on sbottom mass m~ b = m~ bi as in Eq.(6). For tan/5 and 
the mixing parameter M LR) we restrict them to the range tan/5 < 0.25[6], M LR < 
3mg [16]. Other input parameters are vrtz = 91.188GeV, a em = 1/128.8, and Gp = 
1.166372 x \Q- 5 (GeV)- 2 . m w is determined through [17] 

2 / 1 m W\ na ^ /or\ 

"^-Mf^ ^l-Ar ' (36) 
where, to order O (amf / m^) , Ar is given by [18] 

The relative correction to the hadronic cross section as a function of sbottom 
mass is presented in Fig.l for tan/5 = 0.25 and tan/5 = 1, respectively. Since the 
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correction is proportional to 1/ sin 2 /5, the size of correction for tan/5 = 0.25 is much 
larger than the corresponding size for tan/5 = 1. In the range m b < 150 GeV, the 
correction is very sensitive to sbottom mass. The correction can be either negative 
or positive, depending on sbottom mass. For m b > 200 GeV, the correction drops to 
about zero size, showing the decoupling behaviour of MSSM. Each plot in Fig.l has 
a sharp dip, which occurs at the threshold point m t = m b + M y The chargino masses 
M 12 = (230, 100)GeV for tan/5 = 0.25 and M 1>2 = (220, 120)GeV for tan/5 = 1, thus 
the threshold point locates at m b = 76 GeV for tan/5 = 0.25 and m b = 56 GeV for 
tan (3=1. 

Fig. 2 show the dependence of the relative correction to the hadronic cross section 
on the value of tan/5 for sbottom mass m b = 100 GeV and 150 GeV, respectively. 
The correction is very sensitive to tan/5 in the range tan/5 < 1. When tan/5 << 1, 
the correction size gets very large since it is proportional to 1/ sin 2 /5. 

Fig. 3 is the plot of the relative correction to the hadronic cross section versus 
the top-squark mixing parameter M LR . The corresponding neutralino masses in this 
figure are (122, 115, 77, 229) GeV and chargino masses are (230,100) GeV. The 
starting point M^r = correspondes to no-mixing case, at which top-squark masses 
m ti2 = m t L r = (213,216) GeV and the mixing angel 9 = 0. As M LR increases, the 
mass splitting between two top-squarks increases. At M LR = 200 GeV, top-squark 
masses 2 = (103, 285) GeV and the mixing angel 9 = 0.775. The two sharp dips 
in the plot correspond to two threshold points at about M LR = 200 GeV and 230 
GeV, at which m t = + Mqj. 

So, from Figs. 1-3 we found that only for tan/5 < 1 and m b < 150 GeV the 
correction size may exceed 20%. For tan (3 > 1 or m~ b > 150 the correction size can 
only reach a few percent. In Ref.[6], the Yukawa correction from the SUSY Higgs 
sector to the hadronic cross section was found to be very small for tan/5 = 1, on the 
order of a percent, and only for minimum value tan/5 = 0.25 the correction can be 
above 10% but never exceed 20%. Therefore, the genuine supersymmetric electroweak 
corrections are comparable to the Yukawa correction from the SUSY Higgs sector. 

In conclusion, we presented the anlytical expression for the genuine supersym- 
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metric electroweak corrections of order amf/my V to top quark pair production at the 
Tevatron with the consideration of top-squark mixing. Numerical examples showed 
that only for tan f3 < 1 and < 150 GeV the correction can exceed 20%. In the most 
favorable case, these supersymmetric corrections, combined with Yukawa corrections 
of the Higgs sector and also SUSY QCD correction, are potentially observable at 
Tevatron and could be used to place restrictions on MSSM. 
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Figure Captions 



Fig.l The plot of relative correction to the hadronic cross section versus sbottom 
mass, where Mlr = l.bm~ b . 

Fig. 2 The plot of relative correction to the hadronic cross section versus tan (3, 
where Mlr = 1.5m~ b . 

Fig. 3 The plot of the relative correction to the hadronic cross section versus the 
top-squark mixing parameter Mlr. 
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